arXiv: 1501.03157v2 [hep-ph] 26 Mar 2015 


LAPTH-004/15 

MSUHEP-150113 


Lab-frame observables for probing the top-Higgs 

interaction 


Fawzi Boudjema a , Rohini M. Godbole b , Diego GuadagnolF, Kirtimaan A. Mohan c 

“ LAPTh, Universite de Savoie et CNRS, BP110, F-74941 Annecy-le-Vieux Cedex, France 
b CHEP, Indian Institute of Science, Bangalore 560012, India 
c Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, 

USA 

E-mail: fawzi.boudjema@lapth.cnrs.fr, rohini@cts.iisc.ernet.in, diego.guadagnoli@lapth.cnrs.fr, 

kirtimaan@pa.msu.edu 

We investigate methods to explore the CP nature of the tih coupling at the LHC, 
focusing on associated production of the Higgs with a tt pair. We first discuss the 
constraints implied by low-energy observables and by the Higgs-rate information from 
available LHC data, emphasizing that they cannot provide conclusive evidence on the 
nature of this coupling. We then investigate kinematic observables that could probe 
the tih coupling directly, in particular quantities that can be constructed out of just 
lab-frame kinematics. We define one such observable by exploiting the fact that tt spin 
correlations do also carry information about the CP-nature of the tih coupling. Finally, 
we introduce a CP-odd quantity and a related asymmetry, able to probe CP violation 
in the tih coupling and likewise constructed out of lab-frame momenta only. 
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1. Introduction 


The 7-8 TeV runs of the LHC have led to the discovery of a scalar particle with a mass 
m.h — 125 GeV [1-4]. The properties measured so far show very good consistency with those 
expected for the Standard-Model (SM) Higgs boson. Further, these runs have not revealed 
the existence of new particles. The fact remains that the SM cannot address a few pressing 
questions, such as the baryon asymmetry of the universe, the large mass hierarchy in the 
fermion sector as well as an explanation for the Dark Matter abundance in the Universe. 
These issues call for new physics (NP) beyond the SM. Furthermore, the observation of a 
125 GeV elementary scalar, as well as the absence so far of NP at the TeV scale, leave 
unanswered the question of why its mass mh is so different than the gravitational scale. 

In order to ease the explanation of the observed baryon asymmetry, new sources of CP 
violation are desirable. Such sources exist in many simple extensions of the SM. One notable 
example is an extended Higgs sector such as a two Higgs doublet model. Therein, CP 
violation is incorporated in the Higgs sector through mixing of CP-even and -odd states. 
Within these models, the 125-GeV boson identified with the Higgs can have indefinite CP 
quantum numbers due to mixing of CP-even and -odd states. A determination of the CP 
nature of this particle and its interactions may thus hold clue of NP. 

A program to probe the CP nature of the discovered Higgs scalar is already under way 
at the LHC experiments. The pure pseudo-scalar hypothesis has already been ruled out at 
greater than 95% confidence level (CL) and consistency with the CP-even nature established 
by ATLAS [5] and CMS [6,7]. This has been achieved by an analysis of the hZZ coupling 
using (h —y Z^Z^) decay channel. 1 It should be noted however that tree-level coupling 
of the CP-odd component of the Higgs to gauge bosons is in fact not allowed and can only 
proceed through loops. Couplings between the Higgs CP-odd component and gauge bosons 
manifest themselves as operators of dimension six (or higher) in the language of effective 
Lagrangians. The effect of such operators is expected to be suppressed in comparison to 
tree-level interactions. 

On the other hand, the CP-odd component of the Higgs couples to fermions at the tree 
level. As a result the Higgs-fermion couplings provide an unambiguous and more sensitive 
probe of a CP-mixed state compared to Higgs-gauge-boson couplings. 2 It is possible to 
probe Higgs-fermion couplings by studying Higgs decays to fermions. Since these are two- 
body decays of a spin 0 particle, the CP nature of the coupling is reflected in the spin 
correlation of the decay fermions. Luckily, the spin information of the t, r is also reflected 
in the decay products of the same. This as well as their larger couplings, offers possibilities 
of probing the CP nature of the Higgs through an analysis of the hrf and hti coupling. 
Analysis of this coupling using the h —> rf case has been shown to be quite promising for 
this purpose [31-34]. However, a measure of the strength of the coupling CP-odd component 
in the hrr interaction does not automatically qualify as a measure of the same for other 
fermions, i.e. the CP-odd component may not couple to all fermions universally (as is the 
prediction in some NP models). It therefore becomes important to be able to probe the 
CP nature of the Higgs in all its couplings. The largest of all such couplings, tih. cannot 
be tested by direct decay, because h —> tt is not allowed. However, the large value of this 
coupling implies large production rates for associated production of the Higgs with a tt pair, 
and this mode therefore qualifies as the most direct probe of the Higgs-top coupling and of 

1 For discussions of the {h Z 1 '*' 1 Z^*'*) decay mode as a probe of the Higgs CP properties see, for example, 
refs. [8-17]. It is also possible to probe the same in Vector Boson Fusion [18-23] production and associated 
(Vh) production [24-30]. 

2 Note that the Higgs to di-photon decay proceeds through loop processes at LO (unlike decays to W and 
Z bosons), making it sensitive to the CP-odd component of the Higgs. 
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the CP nature of the Higgs. 

Also from a more general perspective, it is well known that the coupling of the Higgs to the 
top quark, is of great relevance, theoretically and experimentally alike. On the theoretical 
side, the importance of the top Yukawa coupling follows from the fact that it is numerically 
very close to unity. Such a large value of the Yukawa coupling is suggestive of an active 
role of the top quark in the generation of the electroweak-symmetry breaking (EWSB) scale. 
As a matter of fact, the Higgs-top interaction has important consequences on spontaneous 
symmetry breaking within the SM - notably, on vacuum stability arguments - as well as 
beyond the SM - where the top drives electroweak symmetry breaking in some scenarios. 
Most importantly, this coupling drives the main production channel at the LHC (gluon 
fusion), and also contributes to the crucial decay of the Higgs into two photons. 

The above considerations justify the importance of measuring the top-Higgs coupling with 
the highest accuracy achievable, and, in particular, of determining it by direct measurement 
via tth production. In this paper we focus on this possibility. 

We parameterize the Higgs couplings to fermions through the effective Lagrangian 

£ hff = -J2-^ h f( a f + ib f^fi W 

f 

where the sum is over all quarks and leptons. In the SM, where the Higgs is a scalar, oj = 1 
and bf = 0 for any fermion /. For a pure pseudo-scalar a/ = 0 and bf / 0. A Higgs 
with mixed CP properties is realized if both aj ^ 0 and / 0. The exact values of these 
coefficients will depend on the specific model. Here we are interested in a model-independent 
approach to determine, from data, the nature of the tth interaction which is potentially the 
largest coupling of all fermions. 

The production and decay rates of the Higgs measured at the LHC [35-37] do provide 
important constraints on the strength of both at,bt . Indirect constraints will be discussed 
in more detail in Sec. 2, with the aim of spelling out the underlying assumptions that enter 
the derivations of these constraints. We will show that strong constraints on at and bt can 
be placed only under these assumptions. 

As argued, the most general and direct determination of the at,bt couplings in eq. (1) 
is possible by measuring tth production. 3 The tth production mode is notoriously hard to 
measure at the LHC, yet feasible. In fact, already with the limited data set of the 7 and 8 
TeV runs of the LHC, the signal strengths in the tth production channel have been measured 
by both ATLAS [48,49] and CMS [50]. Some preliminary studies suggest that a significant 
(> 5a) measurement of Higgs production in the tth channel is possible for upcoming runs of 
the LHC [51-57], 

Needless to say, a measurement of the tth production cross-section alone is not sufficient 
to determine the vertex in eq. (1) completely. To this end, it is necessary to consider in detail 
the tth production and the decay kinematics. In this paper we suggest and discuss useful 
discriminating observables to probe the vertex in eq. (1), with emphasis on those that can 
be defined directly in the lab frame. Note that, on the other hand, we refrain from entering 
the discussion about a precision determination of the vertex. 

The rest of this paper is organized as follows. In sec. 2 we describe and derive indirect 
constraints on the couplings at and b t . In sec. 3 we then proceed to analyse tth production 
at the LHC and construct various observables, including a CP-violating one, that could be 
used to determine the nature of the tth interaction itself. Finally in sec. 4 we summarize 
and conclude. 

3 An alternative approach, which we do not discuss in this work, to study the couplings in eq. (1), is to use 
single top production [38-47]. 
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2. Indirect probes of an anomalous tth coupling 


Electric dipole moments (EDMs) can impose severe constraints on new CP-violating weak 
phases. A scalar with mixed parity that couples to both the electron and the top as described 
by eq. (1) leads to CP violation through interference of the type afbf. Indeed at 2-loop a 
Barr-Zee type diagram induces an EDM for the electron of the form d e oc bt,a e f\ (rnf /ml ) + 
beatf'ii'rn'i /m|), where a e ,b e have been defined in eq. ( 1 ), and fip are known loop functions 
[58]. Under the assumption that the Higgs-electron coupling is standard, a e = 1 ,b e = 0, 
a rather stringent constraint, bt < 0.01, can be realized [59]. Of course, with different 
assumptions on a e ,b e , or even with additional sources of CP violation, this constraint can 
become milder or evaporate altogether. For example, and as emphasized in ref. [59], current 
Higgs data are actually compatible with a Higgs only coupled to third-generation fermions. 
In this case bt values of 0(1) are allowed by the EDM constraints. Furthermore, ref. [60] 
provides another example of multi-Higgs scenario, realized in the framework of a CP violating 
supersymmetric model, in which the current EDM constraints can be satisfied, in spite of 
CP violating couplings between the Higgs states and the top quark. 

It should be noted that, given the smallness of the electron Yukawa coupling, it is unclear 
whether the a e , b e couplings will be accessible experimentally in the near future. In order to 
reconstruct the tth coupling direct probes of the same are necessary, which we will discuss 
in the next section. In this section we focus our attention on the constraints on the tth 
coupling that can be derived from Higgs rate information collected at the LHC. We will 
show that these constraints strongly depend on the nature of the assumption and one cannot 
conclusively determine the tth vertex using signal strengths alone. 


2.1. Constraints from measurements of Higgs rates 

Within the SM, and with Higgs and top masses as measured, there are four main production 
modes of the Higgs at the LHC: gluon fusion (ggF), vector-boson fusion (VBF), Higgs pro¬ 
duction in association with &W/Z boson (VH), and Higgs production in association with 
a tt pair. The gluon-fusion production mode has the largest cross-section at the LHC, and 
the dominant contribution to this process comes from a top loop. The Higgs decay to two 
photons has also a contribution due to a top loop, although the dominant one comes from 
a VU-boson loop. ATLAS and CMS have already put indirect constraints on the value of at 
in eq. ( 1 ), assuming that there are no other sources contributing to the effective couplings 
gg —> h or h —> 77 . At 95% confidence level these constraints read [61,62] 

a t € [-1.2, -0.6] U [0.6,1.3] ATLAS 

at £ [0.6, 1 . 2 ] CMS. 


In this section we extend this analysis by allowing in the fit both at and bt couplings in 
eq. (1), and by including the recently measured tth channel signal strengths [48-50,63,64], 
Higgs couplings to massive gauge bosons are defined by 


C h vv = gm w h UwWnV, + ■ ( 2 ) 

In the SM and at tree level kz = k w = Ky = 1 • As customary, the signal strength measured 
in a particular channel i at the LHC is defined as 


hi = 


n, 


exp 


n 


SM 


(3) 
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where n l exp is the number of events observed in the channel i and ng M is the expected number 
of events as predicted in the SM. In order to contrast specific model predictions with the 
experimentally derived fa we define (as usual) 


Ah = 


n: 


th 


n. 


SM 


^ P a P e p &i 

5> P SM 4 x w ' 


(4) 


Here n l th corresponds to the expected number of events predicted in the hypothesized model 
under consideration; a p corresponds to the cross-section in the p th production mode, i.e. the 
cross-section for Higgs production in one of the four production modes listed earlier; Bt is the 
branching ratio of the Higgs in the ?' th channel; e l p is the efficiency of the p th production mode 
to the selection cuts imposed in the i th channel. Note that the efficiencies in the numerator 
and denominator of eq. (4) are taken to be the same. This is true at leading order for the 
gluon fusion process. 

In order to evaluate the signal strength in the tth production channel, ATLAS [48,49,63,64] 
and CMS [50] first apply some basic selection cuts and then use boosted decision trees (BDT) 
to further separate signal from background. We have checked at parton level that for basic 
selection cuts the efficiency in the two cases of pure scalar vs. pure pseudo-scalar Higgs are 
not significantly different. However, this may not be the case for BDT. We neglect this effect 
here, we namely assume that BDT analyses will have the same efficiency for a scalar and a 
pseudo-scalar Higgs and set them to be equal. 

We next discuss the at and bt coupling contributions to Higgs production from gluon fusion 
and Higgs decay to two photons. The ratio of the Higgs decay width to two photons to the 
SM decay width, at next to leading order and neglecting the small contribution from fermions 
other than the top quark, can be written in the form [65,128] 

T(h -A 77 ) ^ | nwA'fafaw) + (1 - a s /Ti) Affa) | 2 + |§ b t A^(T t )| 2 

T(h -A 77 ) sm I A^(tw) + § (1 - a s /n) A? (r t )| 2 

~ l.df^Kw - 0.21 a t ) 2 + 0.12 bfj . (5) 


Here A* denote the loop functions due to the W loop (A^), the CP-even top coupling (A“) 
and its CP-odd counterpart (A^). The analytic expressions for these functions are given 
in app. B. It should be stressed that, given the measured Higgs and top masses, implying 
Tt = m\/(Am 2 ) « 1 , the top contribution (both scalar and pseudo-scalar) is, to a very 
good approximation, given by its expression in the infinite top-mass limit (see app. B). 
Correspondingly, a s corrections are included in this limit. In the same limit, they affect only 
the scalar contributions, whereas the the pseudo-scalar one is untouched. 

We relate Higgs production through gluon fusion, normalized to the SM value, to the 
corresponding normalized width of Higgs to two gluons. Keeping only the dominant top 
contribution again, we may write , 4 


°{99 h) _ T(h -A gg ) _ 2 . h 2 \ A t( T t)\ 2 (, laA 

°{gg -» h) SM ~ T(h-+gg) SM ~ t t \A?(r t )\ 2 \ 2 vr ) 


~ a 2 t + 2.29 b 2 . 


( 6 ) 


4 The first equality in eq. (6), relating the widths and the cross section, is an exact equality at LO. Luckily, 
in the heavy top-mass limit (m\ 7C 4m[) and because we are considering ratios of cr(gg —> h) and ratios of 
T(ft —¥ gg), the equality holds to a very good approximation also beyond LO. In particular, for the cross 
section the higher-order long-distance corrections involving (infra-red/collinear) emission are universal in this 
limit. There only remains a genuine higher-order correction which depends specifically on the nature of the 
tth coupling. A large part of this finite regular correction cancels when considering the ratios. This explains 
the rather small a s correction that we give in eq. (6). For a thorough analysis see [65,128]. 
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Note that, the indirect effect of the pseudo-scalar contribution in gg — > h (and h — > gg ) is 
more than twice the corresponding scalar contribution (with at = bt). As we will see, in 
direct tth production it is the scalar that contributes the most. 

We now perform a global fit to the Higgs data collected by ATLAS, CMS and Tevatron 
in order to estimate the allowed values of at and bt- We follow closely the procedures of 
refs. [35,36,66]. 

In general, BSM models allow for additional interactions not present in the SM to both 
the scalar and pseudo-scalar components of the Higgs, that may be CP-conserving or not. 
Gluon fusion and Higgs to di-photon decays, being loop-induced processes, are sensitive 
probes of this new physics. In this sense unknown heavy physics not related to the top could 
contribute to the effective operators describing gluon fusion (hG^G^, hG ,Jl, G lw , where G^ 
is the gluon field strength and G^ v its dual) and decays into photons (hF^F^, hF^F^, 
Ffj, u and F IW denoting again the electromagnetic field strength and its dual, respectively). In 
order to account for these additional BSM effects, following ref. [66], we introduce four extra 
parameters K gg , K gg , k 77 and k 77 so that eqs. (5) and (6) are modified as follows 

r 77 /r®“ ^ 1.6^(kw - 0.21 (a t + k 77 )) 2 + 0.12 (b t + k 77 ) 2 ) , 

T gg /r™ ~ (at + K gg f + 2.29 {bt + h gg ) 2 . (7) 

Note namely that the couplings n ggi R gg , ft 77 and k 77 are normalized so that, in these 
observables, they shift at and bt with a relative factor of unity. 

The fit to [it is performed by minimizing the y 2 function defined as 

i X 7 


where are the experimental measurements and <f,; their uncertainties. We take into account 
the possibility of asymmetric errors by using the prescription of ref. [67]. Namely, whenever 
errors are quoted as (Ai)!^, we take at = y if (pH — jli) > 0, and &i = z if (m — fit) < 0 [67]. In 
some of the measured channels, the experimental collaborations have provided information 
on the correlation between different production modes. In this case we modify the y 2 function 
to include these correlations as follows 


X 2 (iJ) 


1 

1 — p 2 






(9) 


where p is the correlation coefficient and i and j correspond to different Higgs production 
modes. The data used in the fits are detailed in app. A. 


2.1.1. Results 

We first perform a fit to the SM couplings a/ and Ky, while setting all other couplings 
to zero. The results of this fit are displayed in fig. 1. Here we only show the contours 
for positive values of ap and Ky- An excess seen initially in the h — > 77 channel (excess 
which is now reduced in ATLAS data and absent in CMS data) pointed to negative values 
of a/, which would have had serious consequences on unitarity [68,69]. In the figure, the 
black dot at (0.97,1.06) indicates the best-fit value, while the yellow, green and blue regions 
correspond to the 68 %, 95% and 99.7% confidence level regions, respectively. The SM value 
of (ny,af) = (1,1) is indicated by a red star. Analyses performed by CMS [70] find a best-fit 
value at slightly smaller values of Ky , while fits performed by ATLAS [71] indicate larger 


6 











K V 

Figure 1: Fit results for a/ vs. Ky. The black dot indicates the best-fit value. The yellow (white), 
green (medium grey) and blue (dark grey) areas represent the 68%, 95% and 99.7% confidence-level 
regions, respectively. The red star shows the SM point («y, o/) = (1,1). 


values of Ky. Since we have used both sets of data, we arrive at a middle point, in very good 
agreement with the SM expectation. We found good agreement with the fits of ATLAS and 
CMS when we use only their respective data sets. 

We next perform a fit to the parameters at and bt - the CP-even and CP-odd Higgs-top 
quark couplings. All other parameters are fixed to their SM values, i.e. Ky = 1, a/ = 1 and 
bf = 0 for any f ^ t. The results of this fit are shown in fig. 2. Similar analyses have also been 
performed in refs. [35,36,72]. We find two best-fit values, (at, bt) = (0.67,0.46) and (at, bt) = 
(0.67,-0.46), which are symmetric about the bt axis as expected. Remarkably, significant 



a t 

Figure 2: Fit results for at vs. bt . Black dots indicate the best-fit values. The yellow (white), green 
(medium grey) and blue (dark grey) areas represent the 68%, 95% and 99.7% confidence-level regions, 
respectively. The red star shows the SM point (a t ,b t ) = (1,0). The fit to a t ,b t is performed while 
keeping all other parameters fixed to their SM values, i.e. Ky = 1, a/ = 1 and bf = 0 for any f ^ t. 
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departures from the SM expectation are still possible for the CP-odd coupling. The shape of 
the 68 %, 95% and 99.7% confidence level regions in fig. 2 can be easily understood by looking 
at eqs. (5) and ( 6 ). In gg -» h the at and bt coefficients enter quadratically, weighed by 
the loop functions A“ and A b t , respectively. Therefore, while gg h production is useful to 
constrain the overall af and b 2 magnitudes, alone it is unable to distinguish between scalar 
and pseudo-scalar effects. Inclusion of the h —> 77 decay channel substantially improves 
the discriminating power. The important point is that, in this decay channel, the scalar¬ 
coupling contribution, contrary to the pseudo-scalar one, interferes with the W contribution. 
In particular, for at > 0, as in the SM, this interference is destructive. On the other hand, for 
Of negative, the branching ratio gets enhanced with respect to the SM one by both the scalar 
and the pseudo-scalar contributions, thus making T(h — > 77 ) too large. This is the reason 
why at < 0 is less favoured than at > 0 in fig. 2. Specifically, a* = 0 does not fit the data 
either because in this case the W loop is too large and cannot obviously be compensated by 
the bt contribution, irrespective of the value of bt■ 

Next we look at the effect of the parameters K gg , k 77 , k gg and k 77 , introduced in eqs. (7). 
In particular, we would like to investigate their impact on the value of bt. By inspection 
of eqs. (7), it is clear that an arbitrary (common) value for k gg and R 77 can always be 
compensated by bt- Therefore, a simultaneous fit of Higgs-rate data to iz gg , k 77 and bt would 
result in the flat direction \bt\ = — \k gg \ = — |/c 77 |, with | 6 t | arbitrary. We quote, as an 
example, a fit where we set R gg = re 77 = — 1, and K gg = k 77 = 0. We find the best-fit point 
(at = 0.67, bt = 1.46) and that the various confidence level contours have shifted by +1 in 
the bt direction, allowing for correspondingly larger values of bt than the fit of fig. 2. We 
would expect a second best-fit point at (at = 0.67, bt = 0.54), according to the discussion 
in the previous paragraph, and as displayed in fig. 2. We actually find that the x 2 value of 
this second solution is not exactly equal to the x 2 a t the best-fit point, although the relative 
difference is puny, 2 x ICC 5 . Exact degeneracy is lifted by the tagged tth data in table 6 , to 
which the ±( 6 f + %g, 77 ) symmetry does not apply, at variance with the rest of the data. This 
example demonstrates that tth data would in principle be able to resolve the degeneracy in 
the bt solutions, but their discriminating power is limited by their small statistical weight as 
compared with the rest of Higgs-rate data. 

Altogether, this example is meant to show the inherent limitation of using indirect effects 
to probe the bt. interaction. As a matter of fact, in spite of using a very minimal set of 
parameters, data does not rule out a non-zero bt . Furthermore, on introducing additional 
sources of pseudo-scalar interactions, even larger values of bt can be accommodated. Finally, 
since signal strengths are CP-even quantities (and therefore not linear in bt), they do not 
provide information on the sign of bt . All such ambiguities in the determination of at and bt 
could only be resolved with more direct probes, as discussed in the remainder of this work. 

3. Associated production of the Higgs with a tt pair 

3.1. Kinematics of tth production: scalar- vs. pseudo-scalar-Higgs cases 

Of the four production modes (ggF, VBF, VH, tth, with V = W ± ,Z) of the Higgs at the 
LHC, tth production has the smallest cross-section. Search strategies for the tth process at 
the LHC have been studied in various Higgs decay modes [73,74]: bb [52,53], t + t~ [54] and 
W + W~ [55-57]. The complicated final state of the process, with the top quark decaying to 
a bottom quark and a W boson, which in turn may decay either hadronically or leptonically, 
as well as the large backgrounds to the process make this a difficult channel to study at the 
LHC. Note, on the other hand, that tth production can be studied very precisely at a future 
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Figure 3: Feynman diagrams for Higgs production in association with a ti pair at the LHC. Diagrams 
where the production is mediated by a Z boson or a photon have been omitted. Three more diagrams 
can be realized by exchanging the two gluon lines in the last row labelled (c). 


linear collider such as the ILC [75]. Sufficiently high rates for this process are possible at such 
colliders [76-82] and can therefore be used to extract CP information [83-91] by exploiting 
angular correlations and/or polarization of the top pair. 

As noted in the previous section, studying tth production at the LHC, though challenging, 
is a necessary undertaking; among the other reasons in order to unambiguously determine the 
parity of the Higgs coupling to the top quark, and to reveal potential CP-violating effects 
in the Higgs-top coupling. In this section we wish to point out the major differences in 
the kinematics of the top and Higgs that a scalar- vs. a pseudo-scalar Higgs entails for tth 
production at the LHC. This has been discussed in the literature in quite some detail. See for 
example refs. [85,92-98] and references therein for studies of the CP nature of the tth vertex 
at the LHC. Many of these employ optimal observables [85,99] or the modern incarnation 
of the technique, the multivariate analysis. The aim of the present work at large, is to 
search for and explore lab-frame observables able to probe the nature of the tth interactions 
at the LHC, in spite of the hadronic environment. Our analyses are performed on 14 TeV 
LHC collisions at the parton level, simulated thanks to the MadGraph package [100]. In view 
of the exploratory nature of this study, we do not consider the effects of NLO corrections, 
backgrounds, hadronization, initial and final state radiation or detector effects. All the 
necessary qualifications will be made as our observables are discussed. 

As a first step, let us try to understand the kinematics of tth production without consid¬ 
ering the decays of the Higgs or the ti quarks. The relevant Feynman diagrams are shown in 
fig. 3. We have grouped the diagrams into three categories: quark-initiated, gluon-initiated 
s-channel and gluon-initiated 4-channel. Diagrams where the production is mediated by a Z 
boson or a photon have been omitted. Three more diagrams can be realized by exchanging 
the two gluon lines in the last row labelled (c). 

The first distribution we consider is the production cross-section near threshold. It has 
been pointed out that the threshold behaviour of the cross-section for a scalar vs. a pseudo¬ 
scalar Higgs is very different at an e + e _ collider [86-88]. More specifically, the rate of 
increase of the cross-section with the centre of mass energy of the collision is suppressed in 
the case of the pseudo-scalar Higgs coupling by a factor of p, where p = (y/s — 2mt — rrih)/y/s 
parametrizes the proximity to the production threshold. This factor can be easily understood 
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Figure 4: (Left panel) The invariant-mass distribution of the tth system, normalized to unity. (Right 
panel) The differential cross-section with respect to the tth invariant mass. In either panel, the SM 
distribution (a* = 1, b t = 0) is shown with a solid black line, the pseudo-scalar case (a* = 0, b t = 1) 
with a blue dashed line, and the CP-violating case (a* = 1,6* = 1) with a dotted red line. 


from arguments of parity and angular-momentum conservation [86]. Close to the energy 
threshold, the simultaneous demand of angular momentum and parity conservation implies 
that for a scalar the total angular momentum of the tth system will be zero, while for a 
pseudo-scalar it will be one. Since the process is mediated through s-channel production, the 
pseudo-scalar production will be suppressed near threshold. Note that the total cross-section 
and not just the behaviour near the threshold is different for a scalar and a pseudo-scalar 
for the same Yukawa coupling strength. 

At the LHC, several competing production mechanisms are at work, and it is non-trivial 
that a similar difference in the threshold rise be also visible. Indeed the same behaviour as 
in the e + e - case is observed in the quark-initiated process of pp collisions, which is a spin-1, 
s-channel process, but this contribution is negligible at the LHC. The dominant gg-initiated 
process, has contributions from both s-channel and t-channel diagrams as shown in fig. 3. 
While for pseudo-scalar production the s-channel displays a similar suppression by p near 
threshold, the t-channel does not. We find however that the cross-section near the production 
threshold in the t-channel displays a suppression by a factor proportional to (m/i/m*) 4 . As 
a result, the production cross-section near threshold does show interesting behaviour. 

In the left panel of fig. 4 we show the normalized invariant mass distributions of the 
tth system for the pseudo-scalar (a* = 0,6* = 1), the scalar (a* = 1,6* = 0) case and the 
CP-violating case. 

We see that the rate of increase of the cross-section with the invariant mass of the tth 
system is much more rapid for the scalar than for the pseudo-scalar case. This is an important 
distinguishing feature and could be used to probe the nature of the Higgs-top quark coupling. 
The right panel of fig. 4 shows the same distributions, but normalized to the total cross- 
section (i.e. da/dM t j h ). We observe, as expected, that for the same coupling magnitude, the 
cross-section for the pseudo-scalar case is suppressed with respect to the scalar case. 

While the invariant mass distribution is a useful observable to probe the nature of the 
Higgs-top couplings, its measurement is not straightforward. In fact, it requires complete 
knowledge of the top and Higgs momenta, whose reconstruction is challenged by uncertainties 
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Figure 5: (Left panel) The distribution of the transverse momentum of the Higgs (py.), normalized 
to unity. (Right panel) The differential cross-section with respect to the transverse momentum of the 
Higgs (pj.). In either panel, the SM distribution (at = 1,6* = 0) is shown with a solid black line, the 
pseudo-scalar case (a t = 0, b t = 1) with a blue dashed line, and the CP violating case (a t = 1,6* = 1) 
with a dotted red line. 


on jet energies and, in particular, by missing energy, in decay channels including neutrinos. 

We note incidentally that, rather than trying to extract the full distribution itself, it might 
be easier to consider ratios of cross-sections in two M**^ intervals. 

The complications mentioned above motivate us to look for alternatives to the invariant- 
mass distribution M t j h . One first possibility, that has also been considered in refs. [92,93], is 
the transverse momentum of the Higgs. Its distributions are shown in fig. 5, with normaliza¬ 
tions analogous to fig. 4. As a general feature, we note that the transverse momentum of the 
Higgs (pt) displays a behaviour akin to the invariant-mass distribution M t j h . Noteworthy 
is the fact that ptf is pushed to larger values in the pseudo-scalar case (a* = 0, 6* = 1) in 
comparison to the SM distribution (a* = 1, 6 * = 0). 

The larger transverse momentum of the Higgs in the pseudo-scalar case will have an effect 
on an observable that can be measured quite easily, namely the azimuthal-angle separation 
between the top quark and anti-quark, A< f>(t, t). In order to measure this quantity one needs 
only to reconstruct one of the top momenta at most. The distribution for this observable 
is shown in fig. 6 for the SM (a* = 1,6* = 0), the pseudo-scalar (a* = 0,6* = 1), and the 
CP-violating case (a* = 1,6* = 1). We see that in either case A(f>(t,t) peaks at large values 
± 7 r. However, for the pseudo-scalar case the distribution is more flat in comparison to the 
SM. This can be understood as follows. For events produced near the energy threshold the 
transverse momentum of the Higgs is small. This means that the top pair will be produced 
mostly back to back. This accounts for the peaks observed at \A(j)(t,t)\ = n. Because the plf 
distribution in the pseudo-scalar case is pushed to larger values, this will give rise to a flatter 
distribution in A (j)(t,t). Considering that the construction of this observable only requires 
information about the direction of the various decay products, it can be readily used in both 
the hadronic as well as semi-leptonic decay modes of the top quarks. Uncertainties in the 
measurement of this observable are likely to be much reduced in comparison to M**^. 

One may also attempt to address the question, which of the observables, M t j h , pi^ or 
A 4>(t,t), better discriminates between scalar and pseudo-scalar production, although at the 
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Figure 6: The distribution of the azimuthal-angle difference between the top pair (A <j)(t,t)), normal¬ 
ized to unity. The SM distribution (a* = 1, bt = 0) is shown with a solid black line, the pseudo-scalar 
case (at =0,6* = 1) with a blue dashed line, and the CP-violating case (a t = 1, b t = 1) with a dotted 
red line. 


experimental level one may rather opt for reconstructing all three and use them in a multi¬ 
variate analysis. To answer this question we perform a likelihood analysis, akin to the 
one described in ref. [28]. For the sake of comparison we assume 100% efficiency in the 
construction of both observables, neglect backgrounds and normalize the total cross-section 
for scalar and pseudo-scalar production to be the same. As a result the luminosities that one 
will achieve from such an analysis are not realistic and are only to be used to appreciate the 
discriminating power of the two observables. We use histograms of the distributions binned 
with 20 intervals in the range (0,2000) GeV, (0,500) GeV and (—7r,7r) for M t i h . pj, and 
A 4>(t,t), respectively. 

In fig. 7, we show the variation of the p-value for the pseudo-scalar hypothesis measured 
from the median value of the SM (null) hypothesis. We have used three likelihood functions, 
L(M t i h ), L(pj) and L(A<p(t,t)). We reiterate that the absolute values of the luminosities in 
this figure are not to be taken seriously, as we are only interested in the slopes of the lines. 
From this figure we can infer that the M t j h distribution has a slightly better discriminating 
power followed by A 4>(t, t) and then by pj,. However, the difference between the three 
likelihoods is very small. Since A cf>(t, t) and p^ will have better reconstruction efficiencies 
and reduced uncertainties in comparison to M t i h . the former are expected to perform much 
better in a more realistic analysis. We conclude that A< fi(t, t) and p j. are better suited 
observables to distinguish between a scalar and a pseudo-scalar hypothesis. 

So far we have only considered the kinematics of tth production, without any regards to 
the decays of the top quarks or the Higgs. Furthermore, the observables we have constructed 
are not directly sensitive to CP-violating effects. We will address these issues in the next 
section. 

It is interesting to note that, although a specific measurement of the tth cross section 
cannot discriminate between a scalar and a pseudo-scalar Higgs, the fact that the distribu¬ 
tions are sensitive to its CP assignment means that by comparing a subset of the same cross 
section one could in principle lift the degeneracy. Normalised to the SM cross section the 
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Figure 7: p- value as a function of the integrated luminosity for the distinction between the SM 
(at = 1 , b t = 0) and the pseudo-scalar case (a t = 0 ,b t = 1 ) using likelihoods constructed from the 
observables M t i h (black solid line), p £ (blue dashed line) and A <j>(t,t) (red dot-dashed line). We only 
include statistical uncertainties associated to the tth signal, in the absence of backgrounds. The two 
horizontal lines indicate the 2a (top line) and 3<r (bottom line) exclusion limits. 


inclusive tth cross section at 14 TeV can be written as 

o-/cr SM ~ a t 2 + 0.42 bf . (10) 

A cut Xcut> such as p 1 ^ > 100 GeV, increases the relative weight of the pseudo-scalar contri¬ 
bution 

^(Xcut)/^ SM (Xcut)/ = a 2 t + 0.60 b 2 t . (11) 

If both these measurements, eq. (10) and eq. (11), were precise enough, combining them could 
return non-zero values for both at and bt ■ While none of these cross sections is a measure of 
CP violation, the combination of both cross sections may lead non zero values for both at 
and bt, which is an indirect measure of CP violation. At the lower centre of mass energy of 
8 TeV, the inclusive cross section benefits less the pseudo-scalar contribution. In fact, even 
at 14 TeV, the pseudo-scalar contribution is enhanced relative to the scalar contribution in 
the more energetic regions of phase space. The cross section at 8 TeV centre of mass is 
parametrised as 

O-STev/o-fTeV ~ a t + °' 31 b t • ( 12 ) 

It should be remembered at this point that precision of cross-section ratios as probes of BSM 
physics is to some extent limited by QCD uncertainties in the cross-section predictions [101]. 

3.2. Spin correlations in tt decay products 

The nature of the Higgs-top coupling in eq. (1) also affects spin correlations between the 
top and the anti-top quarks. The latter can be tested, for example, through azimuthal-angle 
differences between the momenta of the particles involved in the process [102-104], We show 
in fig. 8 the normalized distributions of Acf>(t , t) in unpolarized production for two helicity 
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Figure 8: Distributions of the azimuthal-angle difference between the top pair, A<f>(t,t), normalized 
to unity. The four histograms refer to tth production with like-helicity top pairs (trJ-L + t R t R ) and a 
scalar (black solid line) or a pseudo-scalar (blue dash-dotted line) Higgs, as well as to tth production 
with unlike-helicity top pairs (ti,t R + t R t^) and a scalar (black dashed line) or a pseudo-scalar (blue 
dotted line) Higgs, in the lab frame. 


combinations of the final-state top quarks produced in association with a scalar or a pseudo¬ 
scalar Higgs. The two helicity combinations we consider are like-helicity ( t^L + tutu) and 
unlike-helicity {t R t R + t R t R ) top pairs, in the lab frame. The conventions for helicity states 
and spinors are the same as in ref. [105]. The figure shows that the scalar and especially 
the pseudo-scalar cases produce different effects for different helicity combinations. The 
most striking difference occurs between the unlike-helicity combination for pseudo-scalar 
production, which yields a flat distribution, and the remaining distributions, all clearly 
peaked at \A(f>(t,t)\ = ir. 

A measure of the spin correlations can be defined through the following spin-correlation 
asymmetry in the lab frame 


a{pp ->■ t L t L h ) + a(pp t R t R h) - a(pp t L t R h) - a(pp ->• t R t L h) 
a(pp ->• t L t L h) + a(pp ->• t R t R h) + o(pp t L t R h ) + a(pp ->• t R t L h)' 


(13) 


We find the following numerical values for the spin-correlation asymmetry for the different 
parity admixtures: Ciab(<h = l,b t = 0) = 0.22, Ciab(at = 0,b t = 1) = 0.46 and Ciab(at = 
1, bt = 1) = 0.29. These results can be combined in the following parametric formula 


0.22 of + 0.19 
Clab ~ a 2 t + 0.42 b 2 t 


(14) 


valid for the case of the LHC at 14 TeV. The at,bt dependence of Clab confirms our initial 
remark on the nature of the Higgs coupling affecting spin correlations. To be noted are 
the following points: (i) among the cases considered, the SM predicts the smallest value 
for Clabj (w) although the CP-violating case has a larger value for this coefficient, it is only 
marginally higher than the SM. This is due to the scalar cross-sections being larger than the 
pseudo-scalar ones; (in) the asymmetry in eq. (13) is not sensitive to CP-violating effects as 
it is a CP-even quantity. The same is also true for the observables described in the previous 
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Figure 9: Normalized distributions for A<j> tt (£ + ,£~) in tth production. Distributions are shown for 
the SM (a* = 1 ,b t = 0) (black solid line), for the pseudo-scalar case (a t = 0 ,b t = 1) (blue dashed 
line) and for two CP-violating cases, (a t = 1 ,b t = 1) (red dotted line) and (a t = 1 ,b t = —1) (green 
dot-dashed line). 


section. Note however that a measurement of at and bt is nonetheless an indirect measure of 
CP violation. Theoretically a value for £ which deviates from 0.22 or 0.46, corresponds to 
both at and bt being non zero. 

While the spin-correlation asymmetry in eq. (13) may serve as a yardstick for the order 
of magnitude of the effects to be expected, it is not an easily measurable quantity at the 
LHC. Spin-correlation observables typically exploit the fact that the tt spin information 
is passed on to the kinematic distributions of the decay products of the top quarks. In 
addition, the kinematics of the decay products are more likely to be affected by CP violation 
in the production process than the kinematics of the top quarks themselves, i.e. observables 
constructed using the decay products are more likely to be linearly sensitive to bt-° 

Let us first consider the di-leptonic decay mode 5 6 of the top pair.' It is well known that the 
azimuthal-angle difference between the anti-lepton and the lepton from the decay of t and t, 
respectively, provides a good probe of spin-correlation effects in tt production [104,106-109], 
even in the lab frame. Furthermore, as the lepton angular distribution in the decay of the 
top is not affected by any non-standard effects in the decay vertex, it is a pure probe of 
physics associated with the production process [104,110]. For tth production, the t and t 
are not produced back to back (in the xy plane) since the Higgs momentum adds an extra 
degree of freedom to the system. As a result spin-correlation effects in the azimuthal-angle 
difference will be washed out. It is possible to consider the angles between the two leptons in 
a different reference frame, where the kinematics of the tth system does not dissolve the effect 
of spin correlations. Distributions for such observables can be found in [40,92,111-113]. In 
fig. 9 we show the distribution of one such angle, A cffi(l + ,) [92,111,113]. Aft* (£ + , £~) is 

5 In fact, CP-violating interference terms are more likely to be generated in the matrix element squared 
when we sum over helicities of the decay products since the matrix elements for production and decay can be 
linked through a density matrix. 

6 The observables that we will consider can be altered in an obvious way so that they can be used in the 
semi-leptonic or even hadronic decays. 

' For the di-leptonic channel we apply the following set of cuts: pr of jets > 20 GeV , |rj| of jets < 5 , |ry| of 
b jets < 2.5; pr of leptons > 10 GeV, \r/\ of leptons < 2.5. 
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Figure 10: Normalized distributions for A 9 eh (£ + ,£~) in tth production. Distributions are shown for 
SM (a* = l,6 t = 0) (black solid line), for a pseudo-scalar (a t = 0, b t = 1) (blue dashed line) and for 
two CP violating cases, (a t = l,b t = 1) (red dotted line) and (a t = 1 ,b t = —1) (green dot-dashed 
line). 


defined as the difference between the azimuthal angle of the £ + momentum in the rest frame 
of the top and the azimuthal angle of the £~ momentum evaluated in the rest frame of the 
anti-top [111, 113]. 8 

From the figure we can see that the SM (a* = 1, bt = 0) distribution peaks at Acf> u (£ + ,£~) = 
0, while the pseudo-scalar (at = 0, bt = 1) case has a minimum at A(j> tt (£ + , £~) = 0. We have 
also considered two CP-violating cases, (at = 1, bt = 1) and (at = 1,6* = —1), which show a 
behaviour qualitatively similar to the SM case. Furthermore, since the distributions for the 
two CP-violating cases appear to be the same, we can conclude that the two observables do 
not depend on 6* linearly and hence do not probe CP violation in the production process in 
a direct manner. 

Although the A(f> tt (£ + , £~) and other observables considered in the literature [40] do man¬ 
age to differentiate between a scalar and a pseudo-scalar, they are extremely difficult to 
construct at the LHC, especially because a full reconstruction of all momenta of the tth sys¬ 
tem is necessary. In addition, the uncertainties in the measurement of the various momenta 
involved will carry over to the uncertainties in the measurement of these observables as we 
transform between different frames of reference. We therefore explore the option of con¬ 
structing lab-frame observables. One such observable is A 8 £h (£ + ,£~), defined as the angle 
between the two lepton momenta projected onto the plane perpendicular to the h direction 


8 In constructing the Hr 1 momenta as described, we keep fixed for all events the choice of the x and y axes, 
and the z axis is chosen, as customary, to lie along the beam direction. While individually the azimuthal 
angles for the t + and l~ momenta do depend on the choice of the x and y axes, their difference, as in A</>, 
does not. A rj) depends only on the choice of the beam axis. In fact, one can construct A <f> from the following 
formula 


cos(A <j> tt (l + ,£ )) 


(z x p/_) ■ (z x p/ + ) 

IF/: UFA I 


(15) 


that shows dependence only on the z direction. In this formula, the superscripts t (t) indicate that the given 
momentum is calculated in the rest frame of the t (t). 
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in the lab frame: 


(16) 


cos(A O eh (£ + ,r)) 


(.Ph xpe+)- (Ph x pi-) 
I Ph XP£+\\Ph. Xft- I 


This dehnition can be understood from the following argument. Recall that, for two-body 
tt production, the azimuthal angle of the two leptons is sensitive to spin correlation effects. 
The tth system follows three body kinematics, hence the tt can be understood to ‘recoil’ 
off the Higgs. It follows that, when we project the two lepton momenta onto the plane 
perpendicular to the Higgs direction, the angle between them will also be sensitive to such 
spin-correlation effects. 

The distribution for A 9 ih (£ + ,£~) is shown in fig. 10. From this plot we see that, similarly 
as for the angles considered before, there is an extremum at A 0 ih (£ + ,£~) = 0 for all cases 
considered. The SM distribution displays a pronounced peak at A 9 eh (£ + ,£~) = 0, while the 
pseudo-scalar distribution is smaller and flatter in the whole region [— tt/2, +7t/ 2], whereas 
it is larger at | A9 ih (£ + , £~)\ = -jt. Hence this observable can be used to probe the CP nature 
of the tth interaction. On the other hand, being by its dehnition a CP-even observable, 
A9 eh (£ + ,£~) does not distinguish between the two CP-violating cases (a* = 1,6* = 1) and 
(at = 1,6* = —1), that in fact have exactly the same behaviour in fig. 10. In this respect, it 
is worth noting explicitly that, while the plot in fig. 10 spans the range [—7 r, tt], A 9 eh (£ + ,£~) 
is, according to eq. (16), defined only in the interval [0, n]. In order to assign a given event to 
the [0,7r] or to the [—7r, 0] interval, one needs an observable proportional to sin A 9 ih (£ + ,£~), 
for example sgn (p^ ■ (pg+ x Pe-)), where ‘sgn’ indicates that we consider the sign of the term 
in brackets. 

We conclude this section by noting that, albeit not explicitly shown here, other distribu¬ 
tions that, like the one in fig. 10, are also able to distinguish the different vertex structures, 
would arise if we were to replace one or both of the lepton momenta in eq. (16) by W-boson 
momenta. Such distributions are useful in semi-leptonic or fully hadronic decays of the top 
pair. 


3.3. CP-violating observables 

So far we have confined ourselves to observables that are not sensitive to CP-violating ef¬ 
fects. An observable sensitive to CP violation must be odd under CP transformations. Such 
quantities have been considered in the context of e + e _ colliders [84,88,104,110], and these 
results were exploited in the optimal-observable analysis of ref. [85]. More recently, in the 
context of the LHC, a CP-odd observable was proposed in ref. [40] as follows 

a = sgn [p t tt • (p/J x p/+*)) . (17) 

Here the superscripts indicate that the corresponding momenta are constructed in the centre- 
of-mass frame of the tt system. Because of ‘sgn’, a can only take values of ±1. 

Although this observable is sensitive to CP violation linear in 6*, it suffers from the same 
problem as before: it is very difficult to reconstruct at the LHC as all momenta of the tth 
system need to be determined. We suggest an alternative CP-odd observable that can be 
constructed entirely out of lab-frame quantities: 

P = sgn ((p b - Pi) ■ (pe- x p l+ )) . (18) 

Note that, in order to correctly identify jets originating from a 6 and b quark, one needs 
not reconstruct the top or anti-top momenta, of course. Various algorithms can be used to 
differentiate 6- from 6-jets. 9 

9 See for example refs. [114-117] and references therein. 
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Figure 11: Normalized distributions for (3 ■ A d eh (£ + ,£~) in tth production. Distributions are shown 
for the SM ( a t = 1, b t = 0) (black solid line), for the pseudo-scalar case (a t = 0 ,b t = 1) (blue dashed 
line) and for two CP-violating cases, ( a t = 1 ,b t = 1) (red dotted line) and (a t = 1 ,b t = —1) (green 
dot-dashed line). 


The distribution obtained when we multiply f3 by A 9 ih (£ + ,£~) is shown in fig. 11. This 
distribution displays an asymmetry for the two CP-violating cases. Specifically, the distri¬ 
bution for (at = 1, bt = 1) yields larger values in the positive x-axis, whereas the distribution 
for (at = 1 ,bt = —1) is larger on the negative x-axis. 

We thus have a quantity that not only is sensitive to CP violation but is constructed 
entirely out of lab-frame kinematics. In addition, a measurement of this observable demands 
only reconstruction of the Higgs momentum, whereas reconstruction of the top pair momenta 
is not necessary. Note on the other hand that this observable cannot be generalized easily to 
the case of semi-leptonic or hadronic decays of the top since it is not possible to differentiate 
between the quark and anti-quark jet originating from VF-boson decays. 

It is useful to define CP asymmetries with the observables a x A 9 u (£ + ,£~) [40] and 
0 x A0 ih (£~ ,£ + ) as follows 


a(a X A 9 tI (£+,£~) > 0) - a(a x A9 tl (£+,£~) < 0) 
a (a X A 0 a (£+,£~) > 0) + a(a x A0 a (£+,£~) < 0) 


and 


_ a(0 x A 0 m (r,e+) > 0) -<7(0 x A 9 eh (£-,£ + ) < 0) 
lab “ <r(0 x A0 m (£-,£+) > 0) + <r(0 x A0 eh (£~J+) < 0)' 


( 20 ) 


The dependence of these asymmetries on bt (keeping at = 1 fixed) is shown in fig. 12. We 
observe that both asymmetries are sensitive to the sign of bt (and hence linear in bt), being 
negative for negative values of bt and positive for positive values of this parameter. The 
magnitude of the asymmetry A t t is larger than the magnitude of A i a b for a given value of bt. 
However, we emphasize again that A \ a b is constructed out of lab-frame quantities only and 
as such it is expected to be more easily measurable and to have less systematic uncertainties 
than A t t- 


18 










Figure 12: Dependence of the asymmetries A t i (red solid line) and A\ & \, (blue dashed line) on the 
strength of the CP-odd component bt- In this plot at is kept fixed to unity. 


4. Conclusions 

We have considered a Higgs-top Yukawa coupling that allows for a general scalar and pseudo¬ 
scalar Higgs admixture, and have explored the possibility to probe this coupling in a model- 
independent framework. We find that, although constraining, the information provided by 
the Higgs rates, or by low-energy observables such as EDMs, does not suffice to provide 
conclusive evidence about the nature of this coupling. The arguably best way of probing 
this coupling unambiguously is its direct measurement. While certainly challenging in a 
hadronic environment like the LHC’s, a measurement of this coupling would provide crucial 
information on the properties of the scalar coupled to the SM’s heaviest particle, let alone 
the possibility of unveiling CP-violating effects. 

We have investigated some of the possible kinematic observables that could be used to 
discriminate a scalar- from a pseudo-scalar-like coupling at the LHC, focussing on the possi¬ 
bility of quantities constructed out of just lab-frame kinematics. The information about the 
nature of the coupling is carried by the threshold behaviour of the total invariant mass of the 
tth system, which is however very difficult to reconstruct. We find that similar information is 
encoded in the distributions of two experimentally simpler quantities, namely the transverse 
momentum of the Higgs and the azimuthal-angle separation between the tt pair. 

We furthermore exploit the fact that the information about the nature of the tth interaction 
is also passed on to the decay products of the tt pair. Spin correlations between the t and 
the t are likewise affected by the scalar vs. pseudo-scalar nature of this interaction. We 
suggest several lab-frame observables that are affected by the spin correlations and hence 
can be used to probe the Higgs-top interactions in all possible decay modes of the tt pair: 
di-leptonic, semi-leptonic and hadronic. 

Finally, in the di-lepton channel we construct an observable that bears linear dependence 
on bt and hence is sensitive to CP-violating effects. We determine the corresponding CP 
asymmetry and show how it is sensitive to both the strength and the sign of bt■ 

It goes without saying that, being an exploratory study aimed at the definition of lab-frame 
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observables, the analysis performed here is simplistic. In particular it is a leading-order and 
parton-level analysis. While refinements towards a more realistic analysis (like inclusion of 
NLO, detector smearing, hadronization effects, etc.) will change quantitatively several of 
our distributions, they are not expected to modify our main conclusions. More detailed 
investigations are in progress. 
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A. Data used in fits 

We present in the following tables 1 to 6 the data used in the fits. 


Channel 

Signal strength ji 

m h { GeV) 

Production mode 




ggF VBF VH 

ttH 

ATLAS (4.5fb i at 7TeV + 20.3 ft ” 1 at 8 TeV) [118] 

Inclusive 

1.17 ±0.23 

125.4 

87.5% 7.1 % 4.9% 

0.5% 

CMS (5.1fb -i at 7TeV + lOTfb ” 1 at 8 TeV) [119] 

Inclusive 

1 i4+ u - 2(:i 

124.7 

87.5% 7.1% 4.9% 

0.5% 

Tevatron (lO.Ofb 1 at 1.96TeV) [120] 

Combined 

6-14+3U9 

125 

78% 5% 17% 

- 


Table 1 : Data on signal strengths of h —» 77 recorded by ATLAS and CMS, and at the Tevatron 
The percentages of each production mode in each data are given. 


Channel 

Signal strength /i 

m h ( GeV) 

Production mode 




ggF VBF VH ttH 

ATLAS (4.8fb i at 7TeV + 20.7fb ~ 1 at 8 TeV) [121,122] 

Inclusive 

1 66 +u ' 4t> 
i - DD -0.38 

124.51 

87.5% 

7.1% 4.9% 0.5% 

CMS (5.1fb i at 7TeV + 19.6 fb ” 1 at 8 TeV) [ 6 ] 

Inclusive 

n qo+u. 2 y 
u ‘ yo —0.25 

125.6 

87.5% 

7.1% 4.9% 0.5% 


Table 2 : Data on signal strengths of h — > recorded by ATLAS and CMS. The percentages 

of each production mode in each data are given. 
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Channel 

Signal strength n 

m h ( GeV) 

Production mode 




ggF 

VBF 

VH 

ttH 

ATLAS (25fb 1 integrated luminosity at 7 and 8 TeV. 

) [123] 


ggF 

1 ni+ u - 27 

1 * ul -0.25 

125.36 

100.0% 

- 

- 

- 

YBF 

1 OO+0.53 

1-^0—o -45 

125.36 

- 

100% 

- 

- 

CMS (up to 4.9 fb -i at 7TeV + 19.4 fb” 1 at 8TeV) [15] 

0/1 jet 

0-74™ 

125.6 

97% 

3% 

- 

- 

VBF tag 

0-60±8;^ 

125.6 

17% 

83% 

- 

- 

VH tag 

0.39111? 

125.6 

- 

- 

100% 

- 

WH tag 

0-561^? 

125.6 

- 

- 

100% 

- 

Tevatron (lO.Ofb 1 at 1.96TeV) [120] 

Combined 

0.851]];- 

125 

78% 

5% 

17% 

- 


Table 3: Data on signal strengths of h —> recorded by ATLAS, CMS and Tevatron. The 

percentages of each production mode in each data are given. 


Channel 

Signal strength /i 

m h { GeV) 

Production mode 




ggF VBF VH 

ttH 

ATLAS (4.7ft}- 1 at 7TeV + 20.3ftY J 

at 8TeV) [124] 

VH tag 0.52 ± 0.4 

125.36 

100% 

- 

CMS (up to 5.1fb 1 at 7TeV ± 18.9fb 1 at 8TeV) [125] 

VH tag 1.0 ± 0.5 

125.8 - - 100% 

- 

Tevatron (lO.Ofb 1 at 1.96TeV) [120] 

VH tag 

1 56+ 0 ' 72 

i.OO o 73 

125 

100% 

- 


Table 4: Data on signal strengths of h — ► 66 recorded by ATLAS, CMS and Tevatron. The 
percentages of each production mode in each data are given. 


Channel 

Signal strength /r 

m^GeV) 

Production mode 




ggF VBF 

VH 

ttH 

AT] 

LAS (4.5ft}- 1 at 7TeV + 20.3fb i at 8TeV) [126 


g{ggF) 

H(VBF ± VH) 

L93+];f? 

1 24 +0,57 
1 - Z4 -0.53 

125.36 

125.36 

100% 

59.4% 

40.6% 

“ 

CMS (up to 4.9fb 1 at 7 r . 

reV + 19.7 

Eb i at 8TeV) [127] 


0 jet 

0.34 ± 1.09 

125 

96.9% 

1.0% 

2.1% 

- 

1 jet 

1.07 ±0.46 

125 

75.7% 

14.0% 

10.3% 

- 

VBF tag 

0.95 ±0.41 

125 

19.6% 

80.4% 

- 

- 

VH tag 

-0.33 ± 1.02 

125 

- 

- 

100% 

- 


Table 5: Data on signal strengths of h — > t + t recorded by ATLAS and CMS. The percentages of 
each production mode in each data are given. For ATLAS data we use a correlation of p = —0.5. 
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Channel 

Signal strength [i 

m^(GeV) 

Production mode 




ggF 

VBF 

VH 

ttH 

AT] 

LAS (4.5fb 1 at 7TeV + 20.3fb' 

1 at 8 TeV) [48] 


77 

-i o 2 .t )2 

125.4 

- 

- 

- 

100 % 

ATLAS ( 20.3fb~ i at 8 TeV) [63] 

bb 


125 

- 

- 

- 

100 % 

CMS (up to 5.1fb i at 7TeV + 19.7 fb i at 8 TeV) 

[50] 

77 

o 72.6 

125.6 

- 

- 

- 

100 % 

bb 

+0-7111 

125.6 

- 

- 

- 

100 % 

ThTh 

1 q+6.3 

125.6 

- 

- 

- 

100 % 

4-lepton 

—4 7+ 5 - 0 

^•‘- 1.3 

125.6 

- 

- 

- 

100% 

3-lepton 

+3-1110 

125.6 

- 

- 

- 

100% 

Same-sign 21 

+5-3l?i 

125.6 

- 

- 

- 

100% 


ATLAS ( 20.3 fb" 

L at 8TeV) 

64] 

21epton0r had, 

2^ 

125.0 

- 

- 

- 

100% 

3-lepton 

+2.81U 

125.0 

- 

- 

- 

100% 

21eptonlr/ iarf 

—0 9 +31 
u.y_ 2 0 

125.0 

- 

- 

- 

100% 

4-lepton 

l-Sllo 

125.0 

- 

- 

- 

100% 

llepton2r fta<i 

-9-611? 

125.0 

- 

- 

- 

100 % 


Table 6: Data on signal strengths for various decay modes of the Higgs which is produced through the 
tth production mode, for both ATLAS and CMS. Note that, in the various analyses, contaminations 
from other production modes are negligible. 
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B. Form factors 


The loop functions that we have used in the text appear in the analytic expressions for the 
one-loop induced widths h —> 77 and h —> gg. Keeping only the dominant W and top 
contributions we have (see [65,128]) 


where 


r — 
1 77 — 


r — 
1 gg ~ 


GFa 2 m 3 


128 y/2: 


7 r 


Gpa 2 m. 
64v/27t 3 


H J KW A-w{t W ) + ^ a t A t( r t ) 


+ 


A b (r t ) 


a t A%(r t ) + bt A b t (n) 


( 21 ) 


A a t{r) 

a w( t ) 

A b t {r) 


= -o(t + (t-1 )/(t)), 


= - (2t 2 + 3t + 3(2t-1)/(t)), 

r z 

= 2 /W, 


( 22 ) 


with T* = and 


/( r ) = 


arcsin v'T 


l + \/l—T 1 


i2 


lo Slf7& 


— Z7T 


for r < 1 
for r > 1 


(23) 


Since 77 <C 1, an expansion of the loop functions A^’ b in r* confirms that the departure 
from the infinite mass limit is very small. Indeed we can write 


A t(n) - + ^ + 

A b (r t ) = 2(l + ^ + 0(r t 2 )) . (24) 

With rt = m|/(4mf) ~ 0.15, one sees that corrections to the amplitudes in the 77 —> 0 limit 
are respectively of order 3.5% and 5%. 
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